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Description 

[0001] This invention relates to a method of interpolation in a matrix. The invention is particularly, though not exclu- 
sively, concerned with interpolating between stored values in a memory matrix, for example, for a controller for an 
5 electrical machine. 

[0002] A reluctance machine is a form of electrical machine, consisting of two main parts: the stator, which has salient 
poles with coils around them which typically carry unidirectional current; and the rotor, which has salient poles but no 
windings or permanent magnets. 

[0003] The combined control and power switching electronics, and the reluctance machine are commonly referred 
10 to as a switched reluctance (SR) drive system. The power electronics typically consist of two switches and two diodes 
per phase supplied from a DC source or a rectified AC supply. The switching of these devices has to be synchronized 
to the position of the rotor for the machine to operate correctly. To achieve this, a position sensor, usually having a 
moving part mounted on the rotor shaft of the SR machine, is used to generate an angular rotor position signal which 
is fed back to the control electronics to produce the firing signals for the power devices. A more detailed description 
15 of the switched reluctance machine and its control and power electronics can be found in 'The Characteristics, Design 
and Applications of Switched Reluctance Motors and Drives" by Stephenson and Blake, presented at the PCIM '93 
Conference and Exhibition at Numberg, Germany, June 21-24, 1993. 

[0004] Increasingly, in low-cost, high-volume applications, the control electronics are based around an application 
specific integrated circuit (ASIC) which has a microprocessor embedded in it. 

20 [0005] The control law according to which the SR machine is operated will interpret a specific torque demand on the 
machine to maintain the required machine output. Typically, the controller inputs are the torque demand, or possibly 
a speed demand in a closed-loop speed control system, and the actual machine speed derived from the rotor position 
signal. In response, the controller provides one of a selection of control parameters as outputs for controlling the 
machine. The outputs are commonly an ON command, an OFF command, a freewheel (FW) command and a CHOP 

25 command. While it is possible, in some applications, to programme the controller with the control law in order to compute 
the appropriate output for a given input in real-time, a controller will often employ a look-up table of a set outputs for 
given inputs instead. One form of look-up table comprises a matrix of coordinates in which the inputs of torque and 
speed are used as row and column addresses for accessing an array of stored parameter values. In order to minimise 
the amount of space given over to the matrix, a sparse matrix is sometimes used in which less than a complete set of 

30 outputs is available for a given set of inputs, i.e. there are more discrete combinations of inputs than address locations. 
In this circumstance it is necessary to interpolate between the stored parameter values to achieve smooth control. A 
sparse matrix may be arranged either as a regularly spaced set of parameter values or as an irregularly spaced set of 
parameter values which are clustered around the values of the more common operating conditions. 
[0006] As opposed to determining a parameter value by accessing it in a look-up table, interpolation is a real-time 

35 processing activity. Therefore, there are situations in which the speed of interpolation is a limiting factor on the speed 
of response of the system to an input. It is desirable to be able to reduce the processing time taken by interpolation in 
order to increase the responsiveness of the drive. 

[0007] It is an object of the present invention to provide an efficient method of interpolation between values in a 
matrix in which the amount of processing of data is reduced, 
40 [0008] It is also an object of the invention to provide a method of interpolation between stored values in a processor. 
[0009] Prior art document WO 91/01 527 shows a method of linear interpolation dealing with the problem of speeding 
up the interpolation process by duplicating the hardware. 

[0010] The present invention is defined in the accompanying independent claims. Some preferred features are recited 
in the claims dependent thereon. 

45 [0011] In one form the invention provides a method of interpolation in a processing system between parameter values 
to derive a target parameter value, corresponding to an input ordinate, from an addressable memory device containing 
discrete parameter values at memory locations having at least one dimensional addresses, the method comprising for 
each dimension of address: determining first and second memory location addresses between the contents of which 
is the target parameter value; loading the parameter value held at the first address into a first register: loading the 

50 parameter value held at the second address into a second register; subtracting the contents of the first register from 
the contents of the second register to produce a parameter difference; loading the parameter difference into a third 
register; multiplying the contents of the third register by a gradient value to produce a proportional value; loading the 
proportional value into the second register; and adding the contents of the first register to the contents of the second 
register to produce the target parameter value. 

55 [0012] In another form the invention provides a method of interpolation in a processing system between parameter 
values to derive a target parameter value, corresponding to an input ordinate, from an addressable memory device 
containing discrete parameter values at memory locations having at least one dimensional addresses, the method 
comprising for each dimension of address: determining first and second memory location addresses between the 
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contents of which is the target parameter value; loading the parameter value held at the first address into a first register; 
loading the parameter value held at the second address into a second register; subtracting the contents of the first 
register from the contents of the second register to produce a parameter difference: loading the parameter difference 
into a third register; multiplying the contents of the third register by a gradient value to produce a proportional value; 
loading the proportional value into the first register; and subtracting the contents of the first register from the contents 
of the second register to produce the target parameter value. 

[0013] The invention is particularly applicable to implementation in an ASIC in which the cost of the device is largely 
dependent on the amount of storage and register space it requires. Furthermore, the invention requires only a relatively 
small amount of programme code which reduces storage requirements and increases the speed of response of the 
interpolation. 

[001 4] The invention uses a reduced amount of storage space in the interpolation process by assigning intermediate 
values in the interpolation process to registers to replace accessed stored values no longer required in the calculation 
process. An actual address value for which an interpolated parameter value is required is used to define the parameter 
values in the matrix to either side of the parameter value to be interpolated. 

[0015] In the case of, for example, an electric motor controller, the output from the device could be a control value 
which is dependent on a demand value of speed or torque to the controller and the speed and/or torque feedback from 
the motor itself. In an irregular sparse matrix the regular addresses at which the parameter values are stored correspond 
to irregularly spaced speeaYtorque values. Thus, the speed/torque values represent variable ordinates defining the 
location of the parameter value within the memory device. 

[001 6] The invention can be put into practice in various ways some of which will now be described by way of example 
with reference to the accompanying drawings in which: 



Figure .1 is a schematic block diagram of a switched reluctance drive system incorporating the invention; 

Figure 2 is a block diagram of the main components of the ASIC of Figure 1 ; 
25 Figure 3 is a schematic diagram of some of the components of the ASIC illustrated in Figure 2; 

Figure 4 is a block diagram of the arithmetic logic unit of the ASIC; 

Figure 5 is a graph illustrating chopping current control in a switched reluctance machine; 

Figure 6 is a block diagram of the storage elements of the ASIC of Figure 2; 

Figure 7 is a block diagram of the fetch elements of the ASIC of Figure 2; 
30 Figure 8 is a graphical illustration of a one-dimensional sparse matrix of parameter values addressed according 

to speed values; 

Figure 9 is a graphical illustration of an irregular two-dimensional sparse matrix of parameter values addressable 
according to speed and torque; and 

Figure 1 0 is a graphical illustration of a regular sparse matrix of parameter values addressable according to speed 
35 and torque. 

[0017] Figure 1 illustrates a switched reluctance drive system comprising a switched reluctance motor 10, having a 
rotor 12 and a stator 14. A rotor position transducer (RPT) 16 has an encoded disc 18 mounted on the rotor shaft and 
a fixed single sensor 20. The output of the sensor 20 is arranged to derive rotor position information in accordance 
40 with the movement of the disc 18 as disclosed in EP-A-0630097. A controller 24 has a user demand input 22 in the 
form of either a motor speed and/or a torque demand. The motor 10 is a 3-phase machine, therefore the controller 
generates 6 firing signals for 6 power switches 26. 

[0018] In this embodiment the controller is based on an ASIC and controls the 3-phase motor 10 with the single 
sensor RPT output. Table 1 is a list of the pin descriptions for the ASIC. 



Table 1 : 



50 



55 



ASIC Pin Description 


Pin number 


Description 


1 


VCC1 - Analogue supply +5V 


2 


VCC2 - Digital supply +5V 


3 


RPT - Position sensor signal 


4 


DIR - User direction demand 


5 


SCK - SPI clock and tacho output 


6 


MISO-SPI data in 



3 



EP 0 817 365 B1 



Table 1 : (continued) 



ASIC Pin Description 


Pin number 

r hi ii ui i ii_»sii 


Description 


7 


TEST - Test input 


g 


MOSl - SPI data out 


g 


XTAL1 - Crystal oscillator 


10 


XTAI 2 - Crvstal oscillator 


11 


SS - SPI slave select 


1 c. 


U^RPWM - User interface PWM oulout 


13 


PWMRFF - Chnnninn ctirrent reference 


14 


f^MDI - Dinital suddIv nround 


1S 


f^MDP - Analnnup sunnlv oround 


1 R 


t IQRINJ - 1 kor annlnm >p/Hinital inntit 

UOnllli UOCI Ql IClluUuC'UIUIlCll IIILJUI. 


1 7 


1 RRP.DMP - tntprfarp rnmnaratnr analnnup inont 


1ft 


PDMPRFF - Pnrrpnt nnmnaratnr rpfprpncp inniit 


1Q 


fPOMP - P»-nhasp currpnt comnaralor 




po _, Phasp n hnttnm switch firinn sinnal 


-1 1 


PI - Phn^P P ton Qwitrh firinn sinnal 


22 


BCOMP - B-phase current comparator 


23 


B2 - Phase B bottom switch firing signal 


24 


B1 - Phase B top switch firing signal 


25 


ACOMP - A-phase current comparator 


26 


A2 - Phase A bottom switch firing signal 


27 


A1 - Phase A top switch firing signal 


28 


UVOLT - Under-voltage analogue input 



[001 9] The ASIC may conveniently be divided into three main blocks: the embedded microprocessor 28; the serial* 
peripheral interface (SPI)30; and the main SR control 32. Data can be passed between these blocks by means of 

*o buffers 41 and 43. A block diagram of these elements of the design is shown in Figure 2. This figure also shows the 
main signals which are passed between the blocks. Figure 3 shows the sub-divisions within these main sections. 
[0020] The microprocessor 28 generates the programme counter 33 which is used to address a programme read 
only memory (ROM) 34 within the SPI 30. It also has an output on an address bus 35 which is used to address both 
a data ROM 37 in the SPI, and other storage elements (not shown) within the SR control 32. The other two outputs 

45 from the microprocessor 28 are on a data bus 36 which is buffered in buffer 41 and used as an input 36' to both the 
SPI 30 and the SR control 32, and a set of control signals 38 which, with the address bus 35, determine the destination 
of a given piece of data. Within the microprocessor 28 there are two main sections: an instruction fetch and decode 
section 40 and an arithmetic logic unit (ALU) 42, as shown in Figure 3. The decode section 40 generates the programme 
counter 33 and "fetches" the instruction code and any associated data from the programme ROM 34. It then uses 

so combinational logic to decode the instructions as will be understood by the skilled person. In decoding the instructions, 
control signals are generated and passed to the ALU 42 r which contains four registers A, B. D and M, as shown in 
Figure 4. The registers A, B, D : M are manipulated according to the control signals from the decode section 40. 
[0021] The programme ROM 34 in the SPI 30 is addressed by the programme counter 33, generated by the micro- 
processor 28. The SPI 30 then outputs the program data on bus 44 which is then taken as an input bus 44' to the 

55 microprocessor 28. The data ROM 37 is addressed by a combination of the address bus 35 and the microprocessor- 
generated control signals 38. The data on bus 39 which is extracted from the data ROM 37 is buffered, along with the 
data output 47 of the SR control 32, to form a data input to the microprocessor. 

[0022] There are two sections within the SPI 30: the storage area which contains the programme and data ROMs 
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34 and 37 along with the circuitry required to address them; and input/output (I/O) logic 50 which allows data input to 
the ASIC through to a master-in-slave-out (MISO) pin, and data output through a master-out-slave-in (MOSI) pin (see 
Table 1). These pins may also be addressed by the microprocessor 28 using the address bus 35 and control signals 
38. Data to be output from the ASIC is supplied from the microprocessor data-out bus 36 through buffer 41 to the SPI- 
5 data-in bus 36'. Similarly, data which is input to the ASIC is received by the microprocessor 28 on the microprocessor 
data-in bus 39' via a buffer 43 and the SPI-data-out bus 39. 

[0023] The SR control 32 receives address and control signals 35 and 38, and the data input 36" which is buffered 
from data-out bus 36 of the microprocessor 28. The SR control 32 generates the output on the SR-data-out bus 47, 
which is buffered onto the microprocessor data-in bus 39\ and the firing signals 45 for the six power switches which 
w connect the SR machine to the supply. Figure 3 shows the SR control 32 sub-divided into the following 9 sections: 





1. 


Clocks 70 - 


One system-clock signal is input to the ASIC from a crystal oscillator (not shown). 
Sub-divisions of this clock frequency which are required within the ASIC are 
generated within this block. 


15 


2. 


Frequency Multiplier 72- 


In order to start the motor using the single-sensor coded-disc RPT 1 6, the rotor 
position must be known. To achieve this a "lock" routine is used which puts the 
rotor in a known position when the drive is initially started. To lock the rotor, two 
phases of the machine are energised for a pre-defined time period. One of these 
phases is then switched off so that the rotor is pulled into the position in which it 


20 






is aligned with the energised phase. This is described in EP-A-0630097. The 
circuitry to achieve the timing for this routine is included in this section, as is that 
for a period measurement. The period is read into the microprocessor so that the 
speed of the motor may be calculated for use in the control loop and/or in the 


25 






interpolation routine. 


3. 


RPT decoding 74 - 


The sequence generated by the RPT 1 6 is tested for errors and used to determine 
the active period of each 

phase of the motor as is described in EP-A-0630097. 




4. 


Modes 76 - 


The mode of operation of the motor is determined by the microprocessor, 


30 






according to the speed at which the motor is operating. This will be a choice 








between chopping, single-pulse and continuous current modes, as described in 
EP-A-0534761 . The mode determines the time for which each phase of the motor 
may be active, i.e. the time during which the phase winding is connected to the 
supply. The switching strategy which is to be used is also written from the 


35 






microprocessor 28 and stored in this block. 


5. 


Input-output logic 78 - 


This section includes the storage of values which are written by the 
microprocessor 28 to the SR control 32, and also the buffering of output signals 
onto the SR-data-out bus 47. 




6. 


Angles 80 - 


When the motor is running at relatively high speeds it is operated in the known 


40 






single-pulse control mode. The microprocessor interpolates the data tables to 
determine the angles at which the phases should turn on and off according to 
the speed of and the load on, the motor. These values are 
written to the SR control 32 and stored within the Angles block 80. They are used 
to define the active times of the three phases. 


45 


7. 


Chop 82 - 


The switching strategy which is stored in the Modes block 76 is used to determine 
the firing signals for the power devices in the Chop block 82. The active period 
is either that determined by the RPT decoding block 74 at low speeds or the 
output of the Angles block 80 at higher speeds. The Chop block 82 generates 
the correct switch control signals for the strategy and mode and also implements 


50 






any noise-reduction technique such as ramping down the current at the end of 
the active period according to, e.g. the technique disclosed in US 5461295. 




8. 


Fault 84- 


The firing signals generated in the chop mode are inputs to fault logic. This fault 
block 84 implements drive protection, such as disabling firing when an under- 
voltage condition is detected on the main supply to the drive. 


55 


9. 


Demand Interface 86 - 


The "user" of this drive has a single demand input 22 to control (see Figure 1). 
This input may be an analogue or a digital signal and may 
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(continued) 

be used as a torque, speed or chopping level demand, depending on the control 
strategy which is being implemented. A PWM or analogue demand signal is 
applied to one input of a comparator where it is compared with an appropriate 
reference input signal. The output of the comparator is converted into a n umerical 
value of the error of the demand relative to the reference by an interface unit. 



[0024] Referring to Figure 3 : processing within the ASIC is performed by the microprocessor 28. This microprocessor 
io 28 has the four general-purpose registers A B D and M, contained in the ALU block 42, which may be operated on by 
the instruction set given in Table 2. The ASIC has 8-bit memory addressing. The memory map is defined in Table 3. 
The microprocessor 28 also has a 8-bit index register, contained in block 40. which may be used for indirect loading 
and storing of the general-purpose registers, and an 1 1 -bit start-address register which holds the start-address of the 
main programme code, allowing the use of code which runs only when the device is initialised. 



15 

Table 2: 





Microprocessor Instruction Set 




In^tri if^tinn 

II loll Uvilwl 1 


r)p<5rrintif>n 


Instruction Code (HEX) 


20 I 


ADC 


Arici with Parrv 
rAUvi Willi ually 


40-4F 




ADS 


AHH with ^nturatinn 

rAVJ^J Will) OalUi aUUI I 


60-6F 




PI C 


L/icai octiiy 


1B 




U CO 




F4 


25 


DIV 


Divide 


BO-BF 




HALT 


Qtnn Prnnrammp FypfMitinn 
OLi/kJ n uui aiiii nc i — acouuvi i 


FF 




INC 


Inrrpmpnt thp Inrip* Rpni^tpr 


F3 


30 


JMP 


Jump 


F1 AND F5 


\f\v\r O 




20-2F 




ICR 

Jon 


lurnn tr> ^iihrniitinp 


F? AND F6 




JSRC 


Jump to Subroutine on Condition 


30-3F 


35 


LDD 


Load 


04-0F 




LDDI 


Load Indirect 


9C-9F 




L!M 


Limit 


C4-CF 


40 


MUL 


Multiply 


AO-AF 


NOP 


No Operation 


00 




RTS 


Return from Subroutine 


F0 




SBC 


Subtract with Carry 


50-5F 


45 


SBS 


Subtract with Saturation 


70-7F 




SEC 


Set Carry 


IF 




SHIFT 


Rotate, Arithmetic & Logic Shifts 


80-8B (RIGHT). 90-9B (LEFT) 




SHIFT-N 


Shift by N-bits 


D0-DB 


50 


STR 


Store 


14-17 




STRI 


Store Indirect 


DC-DF 




STRL 


Store on Limit 


10-13 


55 


TJMP 


Jump on Test Condition 


F8 AND FC 


TRD2 


One-dimensional Table Read 


8D 




TRD3 


Two-dimensional Table Read 


8E 
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Table 2: (continued) 




Microprocessor Instruction Set 




Instruction 


Description 


Instruction Code (HEX) 


5 


TSC 


Table Scan 


8C 




TXR 


Transfer Register 


EO-EF 


10 




Table 3: 





ASIC Memory Map 


Address (HEX) 


MEMORY USE 


00-1 F 


MICRO SCRATCH PAD RAM 


30 


INDEX REGISTER 


31 


START ADDRESS REGISTER 


40-4F 


8-BIT SR DATA LATCHED 


50-5F 


16-BIT SR DATA (LATCHED) 


60-FF 


8-BIT CONTROL PARAMETER DATA IN ROM 



Table 4: 



25 



30 



ROM Addressing. 


ADDRESS (HEX) 


DATA 


60-67 


TORQUE 


68-6F 


SPEED 


70-77 


OFF 


78-7F 


FW 


80-BF 


ON 


C0-FF 


CHOP 



40 



50 



55 



[0025] As will be appreciated, the tinning of the power-device switching must be synchronised with the rotor position. 
The switching of the devices is determined by the four parameters, ON angle, OFF angle, free wheel (FW) angle and 
the chopping level. These are shown in Figure 5 of the drawings. They are more fully described in the PCIM'93 article 
referred to above. 

[0026J When the phase switches are ON, both are closed to conduct so that the supply voltage is applied across the 
phase winding. Conversely, the switches are turned OFF when the two switches are opened. A freewheel condition is 
put into effect by closing one switch and opening the other, so that the current in the inductive phase winding recirculates 
and decays more gradually than when both switches are OFF. The chop level is set according to a required output and 
is the upper threshold at which the winding current is maintained in chop control. 

[0027] The control parameters ON. OFF and FW and the chop level are stored in the ASIC in sparse matrices in the 
data ROM 37, along with the torque levels and speeds at which they have been measured. The parameters are ad- 
dressed as 60 to FF hexadecimal (hex), where the addressed matrix ROM contents are as shown in Table 4. Typically, 
the stored values in each matrix are determined empirically by running the motor at pre-determined speeds and loads 
to build up the relevant data. A sparse matrix can be an uneven form of matrix in which the values are clustered around 
the most common operating conditions. Thus, values are often concentrated in the usual operating areas and are 
increasingly sparse in other areas. Nevertheless, in areas of dense, as well as less dense values, interpolation between 
values will be needed to enhance the accuracy with which the system responds to an input. 

[0028] As indicated in Table 4. the values of torque and speed for which measured data have been obtained are 
stored in two contiguous sequences of memory locations which may be thought of as two rows. In this implementation, 
OFF angle and FW angle are taken to vary only with speed and to be independent of torque. This allows the values 
for each of these control parameters to be stored in a single row. By contrast. ON angle and CHOP level are taken to 
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vary with both speed and torque. It follows that storage for these control parameters will require as many rows as there 
are values of torque. 

[0029] In this implementation, eight values of speed and eight values of torque are used and each data value is held 
as an 8-bit number. It will be apparent to the skilled person that other combinations of numbers of data points could 

5 be used, e.g. OFF angle could vary with both speed and torque. 

[0030] Within the SPI 30 there is a sub-level "storage" which contains both the programme ROM 34 and the data 
ROM 37. Figure 6 shows a block diagram of the storage level of the SPI 30. The programme ROM 34 is addressed 
from the microprocessor 28 by Ihe 11 -bit program counter 33. Its output is latched when a programme ROM enable 
signal 92 is high. The data ROM 37 is addressed by the microprocessor 28 via the 8-bit ROM address bus. Although 

10 it is addressed by the microprocessor as 60-FF hex, the data ROM 37 itself has two sets of tables which are addressed 
as 00-9F and A0-FF. A page select system is used to translate the address from the microprocessor 2B into an appro- 
priate form for the data ROM. 

[0031] The microprocessor registers A, B : M and D in Fig. 4 each have three control lines CO, C1 , C2. the functions 
of which are given in Table 5. 

15 



Table 5: 



The Function of the Register Control Lines 


C2 


C1 


CO 


REGISTER ACTION 


0 


0 


0 


DO NOTHING 


0 


0 


1 


LOAD BUS INTO REGISTER 


0 


1 


0 


ROTATE RIGHT 


0 


1 


1 


ROTATE LEFT 


1 


0 


0 


LOGIC SHIFT RIGHT 


1 


0 


1 


LOGIC SHIFT LEFT 


1 


1 


0 


ARITHMETIC SHIFT RIGHT 


1 


1 


1 


ARITHMETIC SHIFT LEFT 



[0032] The registers also have individual SET and clear (CLR) lines as well as carry in (CAIN) and out (CAO) signals. 

The 16-bit outputs of the registers (OUT) are enabled onto the main ALU bus 36 using active-low signals EN AON, 

ENBON, ENDON, ENMON through a 5-input, 16-bit OR gate 98. A 1 6-bit adder 95 is connected to register A through 

a 2-input, 16-bit multiplexer 97. The block FLAG-SET 1 04 sets overflow, zero, and negative flags. 

[0033] Figure 7 is a block diagram of the detail of part of the decode block 40 and shows the means by which the 

instructions are fetched from the program ROM 34 shown in Figure 3. The programme ROM 34 is addressed by the 

programme counter 33. Within the block 40 the instructions are timed using micro-steps which are generated from a 

7-bit m-step counter 110 shown in Figure 7. The counter is reset to 0 before each instruction starts and the active 

signals., during the implementation of the instruction, are timed according to their step location. 

[0034] The programme counter may be loaded using the parallel load (PL) input. The source of the programme 

counter input is a 4-way multiplexer 112, the inputs of which are: 

programme start address 114 - start of the main code loop. 

bus 1 1 6 - if a jump instruction is implemented, the value of the jump destination is written to the programme counter 
via the main data bus. 

return address 118 - if a subroutine is used in the programme then the point to which the programme counter must 
return is loaded via the return address bus. 

fault address 120 - if a fault condition occurs in which the fourth input to the multiplexer is selected, zero is loaded 
into the programme counter. 

[0035] The signals which are used to select the multiplexer inputs, and to parallel-load the programme counter are 
generated by particular instructions, such as the jump and jump-to-subroutine instructions. 

[0036] In order to determine the correct control parameters at a given operating point of torque and speed, the 
microprocessor has to interpolate between the known, stored values. To achieve this, three instructions are designed 
to scan and read the tables, allowing interpolation of the parameters with a minimum of additional code required. 
[0037] The invention is concerned with the interpolation between values in the matrix of stored parameter values. It 
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is desirable in high volume production to be able to use a cost-effective controller. Thus, the ASIC should be designed 
to be as low in cost as possible. Cost savings can be realised by minimising, as far as possible, the amount of data 
storage space required for the programme code and for temporary stores for data manipulation, such as interpolation. 
[0038] The conflicting requirements in a motor controller are for speed of response and the space taken by (and cost 

5 of) the controller components. On the one hand, a look-up table containing a comprehensive range of parameter values 
can be accessed in real-time relatively quickfy as no computation involving the stored parameters need be carried out. 
However to do this requires costly and space-consuming memory. Thus on the other hand, a sparse matrix requiring 
interpolation between stored parameter values requires real-time computation and auxiliary devices for manipulating 
the data. The invention, in one form, seeks to use a sparse matrix while minimising the extent to which auxiliary devices 

10 are required in the manipulation of the data and the processing lime required. This can be done by reducing, where 
possible, the programme code in conjunction with the sparse matrix to reduce the execution time for the interpolation 
routine. According to the invention the requirement for temporary storage in memory and additional registers in the 
execution of the routine can be reduced. 

15 Interpolation of a 1 -dimensional Table 

[0039] As previously explained, the values of the OFF and FW angles are stored in rows in the ASIC as a function 
of speed. Figure 8 shows an example of parameter variation with speed. If the motor is turning at some speed S, the 
table of speed values is scanned to determine the speeds S1 and S2, which are speeds at which values of the parameter 
20 have been obtained and between which the value of S lies. The parameter values P1 and P2 must be extracted from 
the table so that the interpolation may be implemented according to equation (1 ) : 

P=P1 +{[(P2-P1 )*(S-S1 )]/(S2-S1 )} (1 A) 

25 

Where: 

S = the current speed value 
S1 = the lower speed value 
30 S2 = the upper speed value 

P = the interpolated or target value 

P1 = the parameter value at the lower speed value, S1 

P2 = the parameter value at the upper speed value, S2 

35 [0040] It will be apparent to the person of ordinary skill in the art that the interpolation according to equation (1 A) 
can equally well be carried out by starting at the upper speed value and its associated parameter value and progressing 
numerically downwards to the target value. In this case, the equation becomes: 

40 P=P2-{[(P2-Pf)*(S2-S)j/(S2-S1)} (1B) 

The following is based on equation 1A above, but it applies to equation 1B mutatis mutandis. 
Table Scan Instruction (TSC) 

45 

[0041] The table scan instruction TSC in Table 2 is the first part of a torque/speed interpolation routine which imple- 
ments equation (1 ) above. Before the instruction is used, the actual value of speed is loaded into the register A. The 
table scan instruction must have the start address of the part of the table which holds the speed data specified in the 
instruction call. This address is preceded by an e.g. TSC &$68 for the table starting at 68 hex, where $ denotes 
so a hex value. This will be equivalent to TSC &104, i.e. a table starting at 104 decimal. 

[0042] The start address is loaded into an index register located in block 40 to allow indirect addressing of the data 
ROM 37. Successive values of speed are then compared with the value in register A to determine the two speeds in 
the table which lie to either side of the register value. This may be explained in terms of the action on each step: 

55 Step 1 Clear the most significant 8 bits (MSB) of the 1 6-bit bus as the address to be read in is only 8-bits. Enable 
the programme ROM 34 to output the address data. 

Step 2 Enable the address data onto the address bus 35. Parallel load the address data into the index register. 
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Clock a 2-bit counter which is to be used for timing information during step 4. 

Step 3 Put the index register value onto the ROM address bus 35. Again clock enable the 2-bit counter to ensure 
that it is in the correct state at the start of step 4. 

Step 4 Because it is not known how far through the table the desired values will lie : the step counter is disabled, 
during step 4, until either the two values have been located or the end of the table has been reached. 
Therefore., timing during this step is determined from the 2-bit counter which generates three signals, ENO, 
when both outputs are zero, EN1 , when the output is binary 01 , and EN2, when the output is binary 10. It 
is on this last state EN2 (1 0) that the counter is reset, to prevent the fourth possible state 1 1 being reached. 
The three states of the counter are used for each value which is extracted from the table to allow sufficient 
set-up and hold-times for the data ROM addressing. 

At the start of step 4, EN2 is active so the data ROM 37 is enabled, putting the first data value onto the 
bus. The bus value is subtracted from the contents of register A and if the result of this subtraction is 
positive, i.e. if the data value is less than the actual value, the data value is stored, from the bus, into 
register D by setting the control lines of register D to 001 binary, as shown in Table 5. The index register 
is then incremented so that it points to the next data value. The 2-bit counter then resets, returning its 
output to state ENO. During this state, the step counter remains disabled provided that the result of the 
last subtraction was positjve. No other signals are active during the ENO period. When 
EN1 is active, the index register is enabled onto the ROM address bus 35 to set up the address before 
the data ROM 37 is accessed. This pattern continues until the result of a subtraction is negative or the 
end of the table is reached. When either of these conditions occurs, the step counter is re-enabled and 
the instruction moves on to step 5. 

Step 5 Because successive values, less than the required one. have been stored in the register 

D : i.e. when the result of the subtraction was positive during step 4, the lower of the two table points, S1 , 
has already been saved. The value is therefore enabled out of register D and subtracted from the contents 
of register A. This produces the value S-S1 in register A. 

Step 6 The result of the step 5 calculation is enabled out of register A and loaded into register M, ready for the 
division which is to be carried out after the table scan instruction. 

Step 7 The index register was enabled during step 4 to extract from the ROM 37 the value which gave a negative 
subtraction result, i.e. S2. The address to which the index register is still pointing is now enabled onto the 
ROM address bus 35. 

Step 8 The address from the index register is still asserted onto the ROM address bus 35 and the data ROM 37 
is enabled to place the value S2 onto the bus, from which it is taken and loaded into register A. 

Step 9 The index register is decremented to point at the lower of the two table values. 

Step 1 0 The address to which the index register now points is enabled onto the ROM address bus 35. 

Step 1 1 The address from the index register is still asserted onto the ROM address bus 35 and the data ROM 37 
is enabled to place the value S1 onto the bus, from which it is taken and loaded into register D. 

Step 12 The contents of register D is enabled onto the bus 35 and subtracted from the value in register A. The 
result is loaded into register A. 

Step 13 The contents of register A is enabled onto the bus and loaded into register B. Register D is cleared. The 
2-bit counter is reset. The signal, which generates the clock enable for the programme counter and in- 
struction register and resets the step counter, is asserted. 

[0043] The table scan instruction has been designed so that at the end of the instruction, the registers hold the correct 
values to proceed directly with the interpolation calculation. That is, 

B register holds S2 - S1 
M register holds S - S1 
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[0044] Since the evaluation of Equation 1 requires the division of M by B, and since the quotient will be less than 
unity, it is conventional in the art to multiply the dividend (M) by a large constant, say 2 16 befor ethe division is carried 
out. The effect of this constant is taken into account at a later stage in the calculation. 

[0045] In this embodiment, Register D is cleared (i.e. set to all zeros) and appended to register M, effectively providing 
5 multiplication by 2 16 (i.e. 65536). Register B is a 1 6-bit divisor and MD is a 32-bit dividend. The 32-bit quotient is replaced 
in MD and register A is used for the intermediate calculations. An example of the code which might be used in the 
interpolation routine is: 





LDD A #$87 


\ load A with the target value of speed 


10 


TSC &$68 


| table scan from the start address of the speed data 




LDD A &$30 


! toad A from the index register with the address of S1 




SEC 


; set the carry prior to a subtraction - 




SBC #$68 


| subtract the initial table address from the address cf S1 


15 


STR A & TEMPO 


i store the address offset at a temporary location for later use 


DIV U.UL , 


| divide M and D by B to give 65536*(S-S1)/(S2-S1) 




STR D &TEMP1 


; store D, the result of the division, for later use. 



[0046] Table 6 is an example of the table on which the instruction could operate: 

20 

Table 6 



ADDRESS 


68 


69 


6A 


6B 


6C 


6D 


6E 


6F 


DATA (HEX) 


0 


24 


49 


6D 


92 


B6 


DB 


FF 



25 

[0047] Implementing the above instruction steps: 

LDD A #$87 
TSC &$68 

30 index register holds $68, therefore find value 0 from table; 

$87 - $0 = $87 (positive result) therefore store 0 in register D and increment index register to give $69. This retrieves 
value $24 from table; 

$87 - $24 $63 (positive result) therefore store. $24 in register D, increment index register to $6 A & retrieve $49 
from table;. 

35 $87 - $49 = $3E (positive result) therefore store $49 in register D, increment index register to $6B & retrieve $6D 

from table: 

$87 - $6D - $1 A (positive result) therefore store $6D in register D, increment index register to $6C & retrieve $92 
from table; 

$87 - $92 = $FFF5 (negative result) therefore end step 4; 

40 

Calculate S - S1 = $87 - $6D = $1 A and store this result in register M; 

get value S2 (=$92) from ROM and subtract value S1 (stored in register D) from it 

i.e. A holds $92 - $6D = $25 

45 this result is transferred from register A to register B. Register D is cleared (i.e. set to all zeros). 
[0048] The register contents at the end of the instruction are therefore: 

Register A holds 25, the value which is loaded into register B on step 13 
Register B holds 25 = S2-S1 = 92-6D 
so Register M holds 1 A = S-S1 = 87-6D 

Register D holds 0. 

One-dimensional Table Read (TRD2) 

55 [0049] Referring again to Figure 8, the address of P1 is calculated from the address of S1 and pre-loaded into register 
A. This (the address of P1) is enabled onto the bus and parallel loaded into the index register. The address of Pt is 
determined by taking the known address of S1 , subtracting the initial speed table address and adding the resulting 
offset to the start address of the that part of table in which the parameter values are located. 
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[0050] The one-dimensional table read instruction, TRD2, is used, in conjunction with the table scan, when interpo- 
lating data which is dependent on only one variable. Before using the instruction, register A must be loaded with the 
address of P1 . The table read instruction then returns the values of P1 and P2-P1 in registers D and M, respectively, 
for use in the remaining interpolation calculation. The instruction may be explained in terms of the action on each step: 

Step 0 The address value S1 ? which has been pre-loaded into register A : is enabled onto the bus and parallel 
loaded into the index register. 



w 



Step 1 The address held in the index register is enabled onto the ROM address bus. 

Step 2 The address of P1 from the index register is held on the ROM address bus and the corresponding value of 
P1 is enabled out of the data ROM and loaded into register D. The index register is incremented so that it 
holds the address of value P2. 



15 Step 3 Do nothing while allowing set up time for the new address. 

Step 4 The address P2 held in the index register is enabled onto the ROM address bus. 

Step 5 The address from the index register is held on the ROM address bus and the corresponding value of P2 is 

20 enabled out of the data ROM and loaded into register A. 

Step 6 The value held in register D is enabled onto the bus and subtracted from the contents of register A. The 
result, P2-P1, is loaded into register A. 

25 Step 7 The contents of register A are enabled onto the bus and loaded into register M. The signal which generates 
the clock enable forthe programme counter and instruction register and resets the step counter, is asserted. 
The next instruction is enabled out of the programme ROM. 

[0051] An example of the code in which the table read instruction could be used is: 



35 



40 



45 



LDD A #$70 
CLC 

ADC &TEMP2 
TRD2 

STRD&TEMP4 
LDD B &TEMP3 
MUL UL.SL 

CLC 

ADC &TEMP4 
STR A &OFF 



load register A with the start address of the off-angle table 
clear the carry before an addition operation 
add the speed table offset to register A to give the address of P1 
issue table read instruction 

store the contents of register D t P1 , to a temporary location 

load register B with the value 65536*(S-S1)/(S2-S1) generated from a scan of the speed table 
multiply unsigned register B (65536* (S-S1)/(S2-S1)) by signed register M (P2-P1) the result 
being in registers A and D 
clear the carry 

add the stored value P1 to the top 16-bits of (P2-P1)*65536*((S;S1)/(S2-S1) 
store the result of the addition to a location predefined as the off 



Note that discarding the bottom 16-bits of the multiplication result has the effect of dividing by 65536. 



55 



[0052] Table 7 on which this instruction could be implemented is: 

Table 7 



ADDRESS 


70 


71 


72 


73 


74 


75 


76 


77 


DATA (HEX) 


CO 


AO 


AO 


AO 


AO 


9E 


93 


8D 



In this example, the address of P1 is taken to be $74. This is enabled onto the bus in step 0 and loaded into the index 
register. 

[0053] In step 2 the index register is used to address the data ROM. From address $74 we extract the value $A0 
(P1) this is loaded into register D. 

[0054] The index register is then incremented to give the address of P2 (i.e. $75). 
[0055] On step 5 : the value P2 is extracted from the data ROM and loaded into register A. 



12 



EP 0 817 365 B1 



[0056] On step 6 the value P1 (in register D) is subtracted from the value in register A (P2) the result is : 

P2- P1 = $9E-$A0 = $FFFE 

5 

The value is copied into register M on step 7. 

[0057] From the value calculated from the table scan instruction, and the values at the end of the table read instruction t 
the required parameter value (P) can be calculated according to equation (1). The register contents at the end of the 
table read instruction are: 

w 

Register A holds FFFE, the value which is loaded into register M on step 7 
Register B is unused 

Register M holds FFFE = P2-P1 =9E - AO 
Register D holds AO = P1 

15 

[0058] The execution of equation (1 ) then reduces to P = D + M * TEMP3; where TEMP3 holds the result of the table 
scan instruction. 

[0059] It will be appreciated that, while the above descriptions of the Table Scan and Table Read instructions have 
considered that the control parameter was a function of speed alone, the control parameters could equally well have 
20 been a function of torque alone. 

Interpolation of a 2-dimensional Table 

[0060] In this embodiment of the invention, the ON angle and CHOP level are stored as parameters which vary both 
25 with speed and with torque. The appropriate table of data must therefore be scanned to determine the values S1 and 
S2, and T1 and T2, which lie to either side of the measured values of speed, S, and torque, T, respectively. Figure 9 
shows an example of control parameter variation with both torque and speed and the control parameter values, P11 , 
P12. P21 and P22, which must be extracted from the table for the interpolation to be implemented. 
[0061] When the four parameter values have been extracted from the table, the interpolation is performed using 
30 equations (2) to (4) : 

P1=P1 1+[(P12-P1 1)*(S-S1))/(S2-S1) (2) 
P2=P21+[(P22-P21)*(S-S1)J/(S2-S1) (3) 
P=P1+[(P2-P1)*(T-T1)J/(T2-T1) (4) 

40 

where 

P = the interpolated or target value 

P1 = an interpolated parameter value at the lower torque point 
45 P2 ~ an interpolated parameter value at the upper torque point 

P1 1= the parameter value at the lower speed and lower torque point 

P 1 2= the parameter value at the upper speed and tower torque point 

P21 = the parameter value at the lower speed and upper torque point 

P22= the parameter value at the upper speed and upper torque point 
so s~ the actual speed value 

51 = the lower speed value in the table 

52 = the upper speed value in the table 
T = the actual torque value 

T1 = the lower torque value in the table 
55 T2 = the upper torque value in the table 

[0062] Intermediate values can be interpolated from a two-dimensional table in a similar manner according to the 
invention still using the four registers for storing the data and the results of the calculations. 



35 
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[0063] Before the two-dimensional table read can be performed, two table scan instructions are required, one on the 
torque data and one on the speed data. 

[0064] The two-dimensional table read instruction, TRD3 in Table 2, is used, in conjunction with the table scan TSC, 
when interpolating the data. Before using the instruction, register A is loaded with the address of P11 . the table point 
5 which corresponds to the lower values of both torque and speed. The instruction may again be described in terms of 
the action on each step: 

Step 0 Register A, which holds the address of P1 1 , is enabled onto the bus and parallel loaded into the index 

register. 

10 

Step 1 The contents of the index register are enabled onto the ROM address bus. 

Step 2 The address from the index register is held on the ROM address bus and the corresponding value, 

P11 , is enabled out of the data ROM and loaded into register D. The index register is incremented to 
is point to Pt 2. 



Step 3 Do nothing while allowing set up time for the new address. 

Step 4 The contents of the index register is enabled onto the ROM address bus. 

20 

Step 5 The address from the index register is held on the ROM address bus and the corresponding value of 

P12 is enabled out of the data ROM and loaded into register A. 

Step 6 The contents of register D is enabled onto the bus and subtracted from the value held in register A. 

25 The result of the subtraction, Pt 2-P11 is loaded into register A. 

Step 7 The contents of register A are enabled onto the bus and loaded into register M. 

Steps 8-14 Increment the index register seven times to point to P21. (The number seven arises because 8 values 
30 are stored for each control parameter) 

Step 15 The contents of the index register are enabled onto the ROM address bus. 

Step 16 The address from the index register is held on the ROM address bus and the corresponding value, 

35 P21 is enabled out of the data ROM and loaded into register B. The index register is incremented to 

point to P22. 

Step 1 7 Do nothing while allowing set up time for the new address. 

to Step 18 The contents of the index register are enabled onto the ROM address bus. 

Step 19 The address from the index register is held on the ROM address bus and the corresponding value, 

P22 is enabled out of the data ROM and loaded into register A. 

45 Step 20 The contents of register B are enabled out onto the bus and subtracted from the value held in register 

A. The result of the subtraction, P22 - P21 is loaded into register A. The signal which generates the 
clock enable for the programme counter and instruction register and resets the step counter, is assert- 
ed. The next instruction is enabled out of the programme ROM. 



50 

Table 8 



CO 


C1 


C2 


C3 


C4 


C5 


C6 


C7 


C8 


C9 


CA 


CB 


CC 


CD 


CE 


CF 


DO 


D1 


D2 


D3 P11 


D4 P12 


D5 


D6 


D7 
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Table 8 (continued) 



w 



D8 


D9 


DA 


DB P21 


DCp22 


DD 


DE 


DF 


EO 


E1 


E2 


E3 ' 


E4 


E5 


E6 


E7 


E8 


E9 


EA 


EB 


EC 


ED 


EE 


EF 


FO 


F1 


F2 


F3 


F4 


F5 


F6 


F7 


F8 


F9 


FA 


FB 


FC 


FD 


FE 


FF 



Table 8 is an example of an 8-by-8 table starting at address $C0. If P11 is at address $D3 ; increment the index register 
once to get to $D4. Steps 8-14 above increment the index register seven times - it goes from $D4 to $DB : i.e. to point 
at P21. The index register is incremented once more, in step 16, to point to P22, i.e. $DC. 
[0065] An example of the code in which the table read instruction could be used is: 



20 



25 



30 



40 



50 



55 



LDDB&TEMPO 
LSLN #$3 
TXRBA 

CLC 

ADC &TEMP2 
STR A &TEMP2 
CLC 

ADC #$60 
TRD3 

STR D &TEMP4 
STR B &TEMP5 
STR A &TEMP6 
LDD B &TEMP3 



MUL UL,SL 
CLC 

ADC &TEMP4 
STR A&TEMP3 
LDD M &TEMP6 
LDD B&TEMP3 
MUL UL.SL 
CLC 

ADC &TEMP5 
SEC 

SBC &TEMP3 

TXRAB 

LDD M &TEMP1 

MUL SL.UL 
CLC 

ADC &TEMP3 
STR A &CHOP 



load register B with the torque table offset 

shift left by three places which is equivalent to multiplying by 8 

transfer the contents of register B to register A 

clear the carry prior to an addition 

add in the speed table offset to give the position of P11 relative to the start of the data table 
store the relative position of P11 for future use 
clear the carry 

add in the start point of the chop evel data table 
issue table read instruction 

store the contents of register D, Pit, to a temporary location 
store the contents of register B, P21 , to a temporary location 
store the contents of register A s P22-P21 , to a temporary location 
load register B with the value 

65536*(S-S1 )/(S2-S1) generated from a scan of the speed table. This assumes that two table 
scan instructions have already been carried out and the calculations done to generate 65536 
(§i S ih anc * 6 5536 ffp^j the two values then being stored in temporary locations, 
multiply unsigned register B by signed register M(P 1 2-P 1 1 ), the result being in registers A and D 
clear the carry 

add the stored value P11 to the top 16-bits of (P12-P11)*65536*(S-S1)/(S2-S1) 
store the result of the addition, PI, to a temporary location 
load register M with P22-P21 

load register B with the value 65536*(S-S1)/(S2-S1) 

multiply unsigned register B by signed register M, the result being in A and D 

clear the carry 

result of the addition is P2 

set the carry prior to a subtraction instruction 

calculate P2-P1 

transfer the contents of register A to register B 
j load register M with the value 
65536*(T-T1)/(T2-T1) generated from scanning the torque table 
multiply signed register B (P2-P1). by unsigned register M 
clear the carry 

add P1 to (P2-P1)*65536*(T-T1)/(T2-T1) 

store the result of the addition to a location predefined as the chop level store. 



Note that discarding the bottom 16-bits of the multiplication result has the effect of dividing by 65535. 



[0066] The above description concerns interpolation in an irregular sparse matrix in which it is necessary to determine 
the gradient between stored values, i.e. the rate of change between stored values, at given address ordinates, between 
which the interpolation is required. The invention can also be implemented in a regular sparse matrix for which the 
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spacing between the ordinates is constant. 

[0067] Because the spacing between the ordinates is constant, there is no need for the iterative process described 
previously and implemented in the Table Scan Instructions. Instead, the ordinate of, say, speed which lies immediately 
below the specified speed can be identified immediately by simply dividing the specified speed by the spacing between 
the ordinates. For example, Figure 10 shows a regular, 2-dimensionat matrix which has spacing of 16 between the 
points at which the speed is stored. If, say, a value of 71 is presented as the speed value for interpolation of the torque, 
then dividing 71 by 16 gives 4.438. Taking the whole number, 4, indicates that the lower ordinate is the fourth in the 
table above zero, i.e. 4 times 16 = 64. This enables the identification of the lower ordinate to be achieved in a single 
pass through the code for the microprocessor. 

[0068] Using the nomenclature defined above for the irregular matrix, the process for interpolation of a regular matrix 
can be illustrated as follows. 

[0069] A sparse matrix with its addresses is defined as follows: 

msa15, msa14, msa13, msa12 s msa11, msa10, msa09, msa08, P1. PO, S6, S5 : S4, T6, T5, T4; where S6,S5,S4 
are the 3 MSBs of a 7-bit speed signal and T6,T5,T4 are the 3 MSBs of a 7-bit torque signal. P1 and P0 define 
the type of parameters stored in the matrix defined, as before, by: 



P1 


P0 




0 


0 


CHOP 


0 


1 


ON 


1 


0 


FW 


1 


1 


OFF 



msa15...msa08 is the address in the data ROM where the data are located. 
[0070] The equations to solve are : 

Let S = S3,S2,S1 ,S0/((S6,S5 J S4+1)-(S6,S5 J S4)) 
= S3,S2,S1,S0716 



Let T = T3,T2,T1 ,T0/((T6T5T4+1)-(T6,T5,T4)) 
= T3,T2,T1,T0716 



Q=P11 +(P2i-P11)S. 



R=P12 + (P22-P12)S. 



P = Q + (R-Q)T. 



P = P11 + (P21-P11)S+TP12+TS(P22-P12)-TP11+TS(P11-P21) 
= P11 +(P21-P11)S+(P12-P11)T+TS(P22+P11-P12-P21) 

[0071] Before executing the instruction the following data should have been loaded into the 3 registers as shown: 

Register A = = 0,0,0,0,0,0,0,0,0,56,55,54,53,52,51 ,S0 (speed). 
Register B = 0,0,0,0,0,0,0,0,0^6^5^4^3^2^1 ,T0 (torque). 

Register D = msa15,msa14,msa13,msa12,msa11 s msa10,msa09, msaOSP^PO.O.O.O.O.OA 
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Step 0-3 

Logic shift right register A by 4 bits. This discards S3.S2.S1 ,S0. 
Step 4-7 

Logic shift right register B by 4 bits. This discards T3,T2,T1J0. 
5 Step 8-10 

Logic shift left register A by 3 bits. This puts S6,S5 : S4 in the correct location. 
Step 11 

Add register A to Register B. This generates the 6 LSBs of the address. 
Step 12 

10 Add register D to register A. This generates the address of the base point in the table P11 . 

Step 13 

Load the index register with the address of the base point P11 . 
Step 14 

Load the base point of the table, P11 into register D. 
is Step 15-22 

Increment the index register by 8 so that it points at P21 in the matrix. 
Step 23 

Load point P21 into register A. 
Step 24 

20 Subtract register D from register A and put the result (P21-P11) in register A. 

Step 25 

Transfer the contents of register A into register B. (P21 -P11) to B. 
Step 26 

Increment the index register by 1 so that it contains the address of point P22. 
25 Step 27 

Load point P22 into register A. 
Step 28-35 

Decrement the index register by 8 so it contains the address of point P12. 
Step 36 

30 Load point P12 into register M. 

Step 37 

Subtract register M from register A and put the result in register A (P22-P1 2). 
Step 38 

Subtract register M from register A and put the result in register A (P22-P1 2+P11-P21). 
35 Step 39 

Transfer the contents of register A into register M (P22-P12+P11-P21) . 
Step 40 

Load point P12 into register A. 
Step 41 

40 Subtract register D from register A and put the result in register A (P12-P11) 

Step 42 

End of instruction. 

[0072] At the end of the instruction the following partial results have been calculated: 

45 

Register A = P12-P11 
Register B = P21-P11 
Register M = P22-P12+P11-P21 
Register D= P11 

50 

[0073J The result P may now be routinely calculated using existing microprocessor instructions as detailed below: 
Assume S*1 6 is stored at address $S and T*1 6 is stored at address $T. The maximum value of S and T is 1 and thus 
they are stored as 16 times their value. 



STR A & TEMPO j TEMPO contains P12-P11 

STR M & TEMPI ! TEMPI contains P22-P12+P11-P21 

STR D & TEMP2 j TEMP2 contains P11 
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(continued) 





i nn Md9 

L.UU IVI Ot O 


M contains S*16 




MUL.SUUL 


D contains (P21-P11)S*16 


5 


TXR D A 


A contains (P21-P11)S*16 




RORN #$4 


A contains (P21-P11)S 




LDD B&TEMP2 


B contains P11 




L/LO 


Clear carry 


10 








QTD A JtTCIVJlP'? 


1 tZlVl r c. COnidli 15 \r II +■ \n ^ IT II 




1 P\ P\ D9TCRJIDO 

LUU t>& 1 hMrU 


B contains T*16 




LDD M & $T 


M contains P12-P11 


15 


MUL.SLUL 


D contains (P12-P11)T*16 




TXR DA 


A contains (P12-P11)T*16 




RORN #$4 


A contains (P21-P11)T 




LUU d & 1 tlvlr^ 


d coniams r l 1 + [yd. 1 -r 1 i jo 


20 


CLC 






ADC B 


A contains P11 + (P21-P11)S+(P12-P11)T 




STR A&TEMP2 


TEMP2 contains P11+(P21-P11)S+(P12-P11)T 




LDD B &S 


B contains S* 16 


25 


LDD M &T 


M contains T*16 




MUL } UL : UL 


D contains S*T*256 




TXR DA 


A contains S*T*256 


30 


RORN #$4 


A contains S*T*16 


TXR AM 


ivi contains o* i * i o 




LDD B&TEMP 1 


B contains P22-P12+P11-P21 




MUL.SLUL 


D contains {P22-P12+P11-P21)*S*T*16 


35 


TXR DA 


A contains (P22-P12+P11-P21)*S*T*16 




RORN#$4 


A contains (P22-P12+P11-P21)*S*T 




LDD B &TEMP2 


B contains P11 + (P21-P11)S+(Pt2-P11)T 




CLC 




40 


ADC B 


| A contains P1U(P21-P11)S+(P12'P11)T+(P22-P12+P11-P21)ST, i.e. the required value, P. 



[0074] The invention thus provides a method of interpolating between parameter values in a sparse matrix that may 
be regular or irregular. By swapping partial results of the interpolation between registers, minimal use of other storage 
space is needed. The reduction in storage space is accompanied by a reduced programme code. The combination of 

4 5 these two features allows significantly less programme memory to be used and the programme cycle time to be rapid. 
[0075] It will be apparent to the skilled person that the invention can be implemented in processing systems other 
than the ASIC-based SR controller described. Thus, while the present invention is susceptible to various modifications 
and alternative forms, specific embodiments have been shown by way of example in the drawings and description 
herein. However, it should be understood that the invention is not intended to be limited to the particular forms disclosed. 

50 Rather, the invention is intended to cover all modifications, equivalents and alternatives falling within the scope of the 
invention as defined in the appended claims. 

Claims 

55 

1. A method of interpolation in a processing system between parameter values to derive a target parameter value, 
corresponding to an input ordinate, from an addressable memory device containing discrete parameter values at 
memory locations having at least one dimensional addresses, the method comprising for each dimension of ad- 
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dress: 

determining first and second memory location addresses between the contents of which is the target parameter 
value; 

loading the parameter value held at the first address into a first register; 
loading the parameter value held at the second address into a second register; 

subtracting the contents of the first register from the contents of the second register to produce a parameter 
difference; 

loading the parameter difference into a third register; 

multiplying the contents of the third register by a gradient value to produce a proportional value: 
loading the proportional value into the second register; and 

adding the contents of the first register to the contents of the second register to produce the target parameter 
value. 

15 2. A method of interpolation in a processing system between parameter values to derive a target parameter value, 
corresponding to an input ordinate, from an addressable memory device containing discrete parameter values at 
memory locations having at least one dimensional addresses, the method comprising for each dimension of ad- 
dress: 

20 determining first and second memory location addresses between the contents of which is the target parameter 

value; 

loading the parameter value held at the first address into a first register; 
loading the parameter value held at the second address into a second register; 

subtracting the contents of the first register from the contents of the second register to produce a parameter 
25 difference; 

loading the parameter difference into a third register; 

multiplying the contents of the third register by a gradient value to produce a proportional value; 
loading the proportional value into the first register; and 

subtracting the contents of the first register from the contents of the second register to produce the target 
30 parameter value. 

3. A method as claimed in claim 1 or 2 in which the first and second memory location addresses are determined by: 

pointing an index to one memory location in a sequence of memory locations in one part of the memory device 
holding values of ordinates corresponding to the parameter values held in another part of the device; 
loading the input ordinate into the second register; 

loading a first ordinate held at the one memory location pointed to by the index into the first register; 
comparing the ordinates held in the first and second register; 

shifting the index in one direction to point to a next memory location in the sequence if the result of the com- 
parison is not of a predefined character; 

repeating the said comparison of the first and second registers after loading the next ordinate held at the next 
memory location in the sequence into the first register; 

adding a predetermined value to the index pointing to the one memory location address of the first ordinate 
to define the second memory location address in the other part of the device if the result of the comparison is 
of the predefined character; and 

shifting the index in the other direction by one address space to define the first memory address in the other 
part of the device. 

4. A method as claimed in claim 3 in which the gradient value is determined by: 

50 

loading the ordinate at the one memory location in the sequence into the first register; 
loading the ordinate at a previous memory location in the sequence into the second register; 
subtracting the contents of the first register from the input ordinate to produce a first ordinate difference; 
loading the first ordinate difference into the third register; 
55 subtracting the contents of the first register from the contents of the second register to produce a second 

ordinate difference; 

loading the second ordinate difference into a fourth register; 

dividing the contents of the third register by the contents of the fourth register; and 
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storing the result of the division as the gradient value. 

5. A method as claimed in claim 1 or 2 in which the first and second memory location addresses are determined by: 

5 dividing the input ordinate by the value of a spacing between ordinates. corresponding to adjacent parameter 

values, to produce a dividend; 

discarding the least significant bits of the dividend to leave a value comprising the most significant bits; 
shifting an index in one direction by a number of memory location addresses equal to the said value to point 
to the first memory location address: and 
io shifting the index in the one direction by one memory location address to point to the second memory location 

address. 



PatentansprCiche 

15 

1. Verfahren zur Interpolation zwischen Parameterwerten in einem Verarbeitungssystem, um einen Zielparameter- 
wert welchereiner Eingangsordinate entspricht, aus einer AdreBspeichereinrichtung abzuleiten, welche diskrete 
Parameterwerte an Speicherplatzen enthaJL die mindestens eine dimensionale Adresse aufweisen, wobei das 
Verfahren fur jede Dimension der Adresse umfaBt: 

20 

Bestimmen der ersten und zweiten Speicherplatzadressen, zwischen deren Inhalten der Zielparameterwert 

liegt: 

Laden des Parameterwertes, welcher an der ersten Adresse gespeichert ist, in ein erstes Register; 
Laden des Parameterwertes, welcher an der zweiten Adresse gespeichert ist, in ein zweites Register; 
25 Subtrahieren der Inhalte des ersten Registers von den Inhalten des zweiten Registers, um eine Parameter- 

differenz zu erzeugen; 

Laden der Parameterdifferenz in ein drittes Register; 

Multiplizieren des Inhaltes des dritten Registers mit einem Gradientenwert, um einen proportionalen Wert zu 
erzeugen; 

30 Laden des proportionalen Wertes in das zweite Register: und 

Addieren des Inhaltes des ersten Registers zum Inhalt des zweiten Registers, um den Zielparameterwert zu 
erzeugen. 

2. Verfahren zur Interpolation zwischen Parameterwerten in einem Verarbeitungssystem, um einen Zielparameter- 
35 wert, welcher einer Eingangsordinate entspricht, aus einer AdreBspeichereinrichtung abzuleiten, welche diskrete 

Parameterwerte an Speicherplatzen enthait, die mindestens eine dimensionale Adresse aufweisen, wobei das 
Verfahren fiir jede Dimension der Adresse umfaBt: 

Bestimmen der ersten und zweiten Speicherplatzadressen, zwischen deren Inhalten der Zielparameterwert 
40 |jegt; 

Laden des Parameterwertes, welcher an der ersten Adresse gespeichert ist, in ein erstes Register; 
Laden des Parameterwertes, welcher an der zweiten Adresse gespeichert ist, in ein zweites Register; 
Subtrahieren der Inhalte des ersten Registers von den Inhalten des zweiten Registers, um eine Parameter- 
differenz zu erzeugen; 
45 Laden der Parameterdifferenz in ein drittes Register; 

Multiplizieren der Inhalte des dritten Registers mit einem Gradientenwert, um einen proportionalen Wert zu 
erzeugen; 

Laden des proportionalen Wertes in das erste Register; und 

Subtrahieren der Inhalte des ersten Registers von den Inhalten des zweiten Registers, um den Zielparame- 
50 terwert zu erzeugen. 

3. Verfahren nach Anspruch 1 oder 2, bei welchem die ersten und zweiten Speicherplatzadressen bestimmt werden 
durch: 

55 Richten eines Index auf einen Speicherplatz in einer Sequenz von Speicherplatzen in einem Teil der Spei- 

chereinrichtung, welche Werte von Ordinaten speichert, die den Parameterwerten entsprechen, welche in 
einem anderen Teil der Einrichtung gespeichert sind; 
Laden der Eingangsordinate in das zweite Register; 
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Laden der ersten Ordinate, welche an dem ersten Speicherplatz gespeichert ist, auf welchen der Index in dem 
ersten Register gerichtet ist; 

Vergleich der Ordinaten, die im ersten und zweiten Register gespeichert sind; 

Verschieben des index in eine Richtung, urn ihn auf einen nachsten Speicherplatz in der Sequenz zu richten, 
5 wenn das Ergebnis des Vergleiches nicht ein vorbestimmtes Merkmal aufweist; 

Wiederholen des Vergleiches des ersten und zweiten Registers nach dem Laden der nachsten Ordinate, wel- 
che im nachsten Speicherplatz in der Sequenz gespeichert ist, in das erste Register; 
Addieren eines vorbestimmten Wertes zu dem Index, welcher zu der einen Speicherplatzadresse der ersten 
Ordinate weist, urn die zweite Speicherplatzadresse in dem anderen Teil der Einrichtung zu bestimmen, wenn 
10 das Ergebnis des Vergleiches das vorbestimmte Merkma! aufweist; und 

Verschieben des Index um einen AdreGplatz in die andere Richtung, urn die erste Speicheradresse in dem 
anderen Teil der Einrichtung zu bestimmen. 

4. Verfahren nach Anspruch 3, bei welchem der Gradientenwert bestimmt wircj durch: 

15 

Laden der Ordinate an dem einen Speicherplatz in der Sequenz in das erste Register; 

Laden der Ordinate an einem vorhergehenden Speicherplatz in der Sequenz in das zweite Register; 

Subtrahieren der Inhalte des ersten Registers von der Eingangsordinate um eine erste Ordinatendifferenz zu 

erzeugen; 

20 Laden der ersten Ordinatendifferenz in das dritte Register; 

Subtrahieren der Inhalte des ersten Registers von den Inhalten des zweiten Registers, um eine zweite Ordi- 
natendifferenz zu erzeugen; 

Laden der zweiten Ordinatendifferenz in ein viertes Register; 
Teilen der Inhalte des dritten Registers durch die Inhalte des vierten Registers; und 
25 Speichern des Ergebnisses der Division als Gradientenwert. 

5. Verfahren nach Anspruch t oder 2, bei welchem die ersten und zweiten Speicherplatzadressen bestimmt werden 
durch: 

30 Teilen der Eingangsordinate durch den Wert eines Abstandes zwischen den Ordinaten entsprechend den 

benachbarten Parameterwerten, um einen Dividenden zu erzeugen; 

Aussondern der geringstwertigen Bits des Dividenden . um einen Wert zu belassen, welcher die hochstwertigen 
Bits umfaBt; 

Verschieben eines Index in eine Richtung um eine Anzahl von Speicherplatzadressen, die dem genannten 
35 Wert gleich ist, um zur ersten Speicherplatzadresse zu weisen; und Verschieben des Index in die eine Richtung 

um eine Speicherplatzadresse, um zu der zweiten Speicherplatzadresse zu weisen. 



Revendi cations 

40 

1 . Procede ^interpolation dans un systeme de traitement entre des valeurs parametriques pour obtenir une valeur 
parametrique cible, correspondant a une ordonnee d'entree, d'un dispositif de memoire adressable contenant des 
valeurs parametriques discretes a des emplacements de la memoire ayant des adresses au moins unidimension- 
nelles, te procede comprenant pour chaque dimension d'adresse les etapes de : 

45 

determination d'une premiere et d'une deuxieme adresses d'emplacements de memoire entre les contenus 
desquels se trouve la valeur parametrique cible; 

chargement de la valeur parametrique conservee a la premiere adresse dans un premier registre; 
chargement de la valeur parametrique conservee a la deuxieme adresse dans un deuxieme registre; 
50 soustraction du contenu du premier registre du contenu du deuxieme registre pour produire une difference 

parametrique; 

chargement de la difference parametrique dans un troisieme registre; 

multiplication du contenu du troisieme registre par une valeur de gradient pour produire une valeur propor- 
tionnelle: 

55 chargement de la valeur proportionnelle dans le deuxieme registre; et 

addition du contenu du premier registre au contenu du deuxieme registre pour produire la valeur parametrique 
cible. 
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2. Procede ^interpolation dans un systeme de traitement entre des valeurs parametriques pour obtenir une valeur 
parametrique cibie ; correspondant a une ordonnee d'entree, d'un dispositif de memoire adressable contenantdes 
valeurs parametriques discretes a des emplacements de la memoire ayant des adresses au moins unidimension- 
nelles. le procede comprenant pour chaque dimension d'adresse les etapes de : 

5 

determination d'une premiere et d'une deuxieme adresses d'emplacements de memoire entre les contenus 
desquels se trouve la valeur parametrique cible; 

chargement de la valeur parametrique conservee a ia premiere adresse dans un premier registre; 
chargement de la valeur parametrique conservee a ia deuxieme adresse dans un deuxieme registre; 
w soustraction du contenu du premier registre du contenu du deuxieme registre pour produire une difference 

parametrique; 

chargement de la difference parametrique dans un troisieme registre; a multiplier le contenu du troisieme. 
registre par une valeur de gradient pour produire une valeur proportionnelle: 
chargement de la valeur proportionnelle dans le premier registre; et 
is soustraction du contenu du premier registre du contenu du deuxieme registre pour produire la valeur para- 

metrique cible. 

3. Procede selon la revendication 1 ou 2. dans lequel les adresses d'emplacements de la premiere et de la deuxieme 
memoires sont determiners de la maniere suivante : 

20 

on pointe un index sur un emplacement de memoire dans une sequence d'emplacements de memoire d'une 
partie du dispositif de memorisation contenant des valeurs d'ordonnees correspondant aux valeurs parame- 
triques conservees dans une autre partie du dispositif: 
on charge les ordonnees d'entree dans le deuxieme registre; 
25 on charge une premiere ordonnee conservee dans le premier emplacement de memoire pointe par I'index 

dans le premier registre; 

on compare les ordonnees conservees dans le premier et le deuxieme registres; 

on decale I'index dans une direction pour pointer sur un nouvel emplacement de memoire de la sequence si 

le resultat de la comparaison n'a pas un caractere predefini; 
30 on repete ladite comparaison du premier et du deuxieme registres apres chargement de I'ordonnee suivante 

conservee dans ['emplacement de memoire suivant de la sequence dans le premier registre; 

on ajoute une valeur predSterminee a I'index pointant sur la premiere adresse d'emplacement de memoire de 

la premiere ordonnee pour definir la deuxieme adresse d'emplacement de memoire dans I'autre partie du 

dispositif si le resultat de la comparaison a le caractere predefini; et 
35 on decale I'index dans I'autre direction d'un seul espacement d'adresse pour definir la premiere adresse de 

memoire dans I'autre partie du dispositif. 

4. Procede selon la revendication 3, dans lequel la valeur du gradient est determined de la maniere suivante : 

40 on charge I'ordonnee du premier emplacement de memoire de la sequence dans le premier registre; 

on charge I'ordonnee a un emplacement de memoire precedent de la sequence dans le deuxieme registre; 
on soustrait le contenu du premier registre de I'ordonnee d'entree pour produire une premiere difference d'or- 
donnees; 

on charge la premiere difference d'ordonnees dans le troisieme registre: 
45 on soustrait le contenu du premier registre du contenu du deuxieme registre pour produire une deuxieme 

difference d'ordonnees; 

on charge la deuxieme difference d'ordonnees dans un quatrieme registre; 

on divise le contenu du troisieme registre par le contenu du quatrieme registres; et 

on memorise le resultat de la division comme valeur de gradient. 

50 

5. Procede selon la revendication 1 ou 2, dans lequel la premiere et la deuxieme adresses d'emplacements de la 
memoire sont determinees de la maniere suivante : 

on divise I'ordonnee d'entree par la valeur d'un espacement entre les ordonnees correspondant a des valeurs 
55 parametriques adjacentes pour produire un dividende; 

on rejette les bits les moins significatifs du dividende pour laisser une valeur comprenant les bits les plus 
significatifs; 

on decale un index dans une direction d'un nombre d'adresses d'emplacements de la memoire egal a ladite 
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valeur pour pointer sur la premiere adresse d'emplacement de la memoire; et 

on decale I'index dans la meme direction d'une adresse d'emplacement de memoire pour pointer sur la deuxie- 
me adresse d'emplacement de memoire. 
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